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The general problem is considered of calculating the voltage when 
a sinusoidal current generator is connected to a parallel combination 
of a linear and a nonlinear resistance load. A practical algorithm is 
described for computing any desired moment and any Fourier com- 
ponent of voltage. An alternative approximate treatment is also pre- 
sented which avoids numerical integrations and is valid when the 
nonlinear characteristic is rapidly varying. Both methods are applied 
to a photovoltaic optical telephone receiver employing a silicon nirp- 
photodiode and a conventional ring armature telephone receiver cou- 
pled by a transformer. Harmonic distortion is presented for several 
illustrative cases. A clipping level is defined for the receiver, and it is 
proposed that the receiver clipping level should be matched to the 
clipping level of the analog optical channel bringing the signal. On the 
basis of this principle a simple procedure is given, along with the nec- 
essary curves, for determining the required optical power at the source 
and the optimum transformer ratio for any value of transmission loss 
between source and receiver. An illustrative example is given for an 
analog dynamic range of 18 dB that requires a peak source power in the 
lightguide of 0.9 mW. This type of receiver may find limited application 
if lightguides ever serve customers directly. 

I. INTRODUCTION 

There is now strong indication 12 that optical transmission using 
lightguides 3 and optical cables is technologically approaching a readiness 
for use in telecommunications. While we are not here suggesting that 
any extensive application of optical telephones is foreseeable, we have 
nevertheless found the prospect of limited special applications suffi- 
ciently interesting to undertake a study of optical telephone receivers 
from a device point of view. The receiver is only one of a number of de- 
vices, some of which perhaps have not even been invented yet, that would 

1569 



be required in an optical telephone. In addition, further devices would 
be required in the electrical-optical interface between the lightguides 
and the metallic network. The receiver, however, determines one key 
property of the system, the optical power required at the interface to 
transmit speech to the human ear with an acceptable volume and qual- 
ity. 

We presuppose that an optical telephone receiver is required to con- 
vert analog-modulated light power to sound pressure at the ear with no 
other power available. There are two mechanisms that might be em- 
ployed to do this using analog modulation: the optoacoustic effect in 
which sound is directly produced when power-modulated light is ab- 
sorbed, and the photovoltaic effect in which an intermediate electrical 
signal is produced which produces sound by way of an electrical ear- 
phone. We have previously completed a theoretical 4 - 5 and experimental 6 
study of optoacoustic receivers in which it was necessary to solve a variety 
of linear acoustical problems to establish the feasibility of the device and 
to obtain its response. Nonlinear distortion was not considered and is 
not believed to be very important for optoacoustic receivers. Optimi- 
zation in that case is a matter of maximizing the response subject to the 
requirement of a flat response over the telephone voice band, 300- 
3300 Hz. 

Subsequently we have attempted to define an optimum photovoltaic 
receiver in a similar way on the basis of maximizing the linear response. 
It might at first appear that this is a simple problem, because a suitable 
photodiode and earphone are already available, the required frequency 
response has already been engineered into the earphone, and one might 
expect that it is only required to match the impedances of the photodiode 
and the earphone by a simple two-port network (e.g., a transformer). 
However, we have now concluded that no optimum of this type exists 
for the photovoltaic receiver, and a different principle of optimization 
is required which is based on the nonlinearity of the photodiode and the 
quality of speech reproduction required in the system. Thus it has not 
been possible to keep system concepts completely out of the discus- 
sion. 

We propose here a very simple principle of optimization, which we call 
"quality matching," which leads to important conclusions about the 
optical power at the interface, the dynamic range of the system, sound 
levels in the system, and of course, the receiver design. 

We base our discussion on the circuit of Fig. 1. Modulated light power 
u is conducted by a lightguide to the photodiode, which may be regarded 
as a current generator g in parallel with a junction current j(v). The 
photodiode assumed is a specific silicon n + irp + structure 7 designed and 
packaged for lightguide use with a light-sensitive area of diameter 
80 nm and a quantum efficiency at 900 nm of about 77 = 0.8. It is ideal for 
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Fig. 1 — Circuit of photovoltaic receiver. The earphone resistance p is transformed 
tor. 

our purpose because of its low series resistance (« 1 ft). The earphone 
is the ring armature telephone receiver 8 which provides good sensitivity 
with essentially flat response over the band 300-3300 Hz. The photo- 
diode and earphone are coupled by an ideal transformer having a turn 
ratio n adjustable by virtue of primary taps. The size and cost of the 
transformer would be approximately proportional to the maximum value 
of n. The effective load resistance seen by the photodiode is zero at dc 
and r = n 2 p at all signal frequencies, p being the earphone resistance. 
The value of n is to be selected at the time of installation in accordance 
with our optimizing principle. 

The first result to be described in this paper is the nonlinear analysis 
itself, in Section II. Mathematically our problem is the following: Find 
the periodic voltage response v(t) to a source current g\ cos cot in Fig. 1 
assuming r is linear and independent of frequency (except dc), j(u) is 
nonlinear and monotonic, and v{t) has zero average value. What makes 
the problem awkward to treat by textbook methods 9 is the restraint on 
the average value. The analysis, so far as we know, is not covered in texts 
on nonlinear circuits, and may be applicable to a variety of situations. 
In Section III is given an approximate method called the clipping model 
which we have found quite reliable and especially appropriate for the 
photovoltaic receiver. 

The main body of the paper consists of Sections IV, V, VI, and VII, 
devoted to the sensitivity, harmonic distortion, clipping level, and quality 
matching of the receiver, respectively. The sensitivity is an inverse 
measure of response (the smaller the better!) defined here in the same 
way as in our previous work as the amplitude of sinusoidal (power) 
modulation of u required to produce at the ear the average speech power 
level (81 dB SPL) found in the telephone network surveys. It is interesting 
to note that the minimum sensitivity is achieved for r « 3 X 10 5 ft, which 
is considerably smaller than the small signal resistance of the photodiode, 
w 8 X 10 9 ft. This shows the essential importance of the nonlinear analysis 
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presented here since a linear analysis would lead to the equality of these 
two resistances. The minimum sensitivity is not, in general, the optimum 
because of nonlinear distortion. 

We have made extensive calculations of the second, third and fourth 
harmonic distortions and the total harmonic distortion. Some curves 
of total harmonic distortion for a few selected cases are presented in 
Section V. We have found it difficult, however, to draw concrete con- 
clusions from a consideration of the harmonic distortion that could be 
used as the basis for an optimizing principle. Rather, we have turned to 
clipping as the most convenient, relevant, and useful way of specifying 
the nonlinear distortion. 

The clipping model, Section III, is based on the assumption that the 
distortion of the waveform is an abrupt one-sided clipping of the peak. 
This assumption is shown to be an adequate modeling of the nonlinear 
effects of an exponential junction characteristic. It is shown in Section 

VI that a clipping level can be defined which is analogous (except for 
being one-sided) to the clipping level of an analog-modulation channel. 
Our optimizing principle, "quality matching," then follows in Section 

VII in an obvious way, namely that the clipping levels of the receiver and 
the analog light channel feeding the receiver should be set equal to each 
other. The receiver then has the minimum sensitivity consistent with 
the requirement that the quality of the channel not be degraded by the 
receiver. This means that the system could operate at the minimum 
power consistent with a given dynamic range. By using the curves given 
in Section VII, quality matching can be carried out to determine for any 
desired dynamic range the optical power required and the correct 
transformer ratio for the receiver needed for the particular value of the 
transmission loss at that receiver. 

A summary with discussion and conclusions is given in Section VIII. 
This is followed by two appendixes containing useful background ma- 
terial on speech quality and the electrical-optical interface. 

II. NONLINEAR ANALYSIS 

In the circuit of Fig. 1 ignore the light, the transformer, and the ear- 
phone; assume that the effective load resistance r is constant at all signal 
frequencies of interest and zero at dc. Let the current generator be 

8 = So+gi cos 
6 = cot 

g o >gi>0. (1) 

Define averages over time by 

W))6 = ^- C f(6)dd^f . (2) 
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The equation of the circuit is 

v - rj = rgi cos 9 - rj(v), (3) 

which implies the relation 

v = 0. (4) 

Concerning ;'( v) we only assume 

j(uY = dj/dv > 0, y(O) = (5) 

which rules out negative-resistance instabilities. We call v{6) the solution 
and 0(") the inverse solution. 
Define the parameters 

v(ir) = a, u(0) = b, (6) 

and 

q ± ^[(b±a) + r(j(b)±j(a))]. (7) 

We can assume v and 6 are restricted to the domains 

a < v < b, O^d^ir (8) 

and the inverse solution 0{v) is unique. It follows that 



and 



1 r b 
v Q = a + - I 6(p) dp 

■K J a 

jo=j(a) + - C b j(uye( v )d v . 

7T %)a 



0) 



(10) 



From (3) we obtain the inverse solution 

d(v) = cos" 1 [(i' + rj(v) - q+)/q-] (11) 

and the implicit relations which determine a,b 

qjjo = r (12) 

(q-/q + )j = gi. (13) 

Our algorithm proceeds as follows: (i) choose initial values for a,b; 
(ii) compute q± from (7) and; from the integral (10); (Hi) test (12) and 
(13); (iv) iterate this procedure with adjusted values of a,b until the 
desired precision is achieved; (u) the inverse solution is now given by (11) 
with the final values of a,b; (vi) as a check, evaluate v Q from (9) and test 
(4). Since (9) and (10) involve numerical integrations, this algorithm 
requires a large computer. 
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A function of v such as j(v) can be expanded in a Fourier cosine se- 
ries 

MB)) =j + £ j k cos kd (14) 

k=i 

j k = 2U(u(e))coakd) e . (15) 

Write (15) in the form 

J*"T f JW Siak0 dv. (16) 

■KK %)a 

It follows that the Fourier coefficients of v are 

2 r b 
p* = — I sin kd du. (17) 

TTK Ja 

From (3) it follows that Vk is also given by 

v k = q-hi ~ rj k (k = 1,2, . . .), (18) 

which is preferable to (17) for numerical work because it automatically 
becomes exact as jk -* 0. The moments of v{6) can be obtained in the 
same way as (10) 

(v(0) k ) 9 -a k + - C b v k ~ 1 e{v)dv. (19) 

■K J a 

Define a normalized waveform 

0(0) = („(0) - q+)/q- (20) 

which has the property of reducing to the input waveform cos 6 whenever 
j{v) = 0. If j(v) is a rapidly increasing function of v (such as the expo- 
nential characteristic of a junction), it may be that the nonlinearity of 
j(p) can be neglected over part of the cycle while over the rest of the cycle 
the nonlinearity effectively clamps v at a constant upper limit. This is 
one-sided abrupt clipping. It is customary to define the clipping fac- 
tor $ in terms of the ratio of the true analog peak to the clipped peak as 
follows 

* = -2Olog0 max , (21) 

where from (20) 

ma x = 0(0) = (b - q+)/q-. . (22) 

The power delivered to the load at frequency ku is 

Pk = pi/2r. (23) 
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A common way of specifying nonlinear distortion is in terms of the ra- 
tios 

d k =Pk/pi (fe = 2,3,...) (24) 

and 

d= Ed k = (2(u 2 ) g /p 2 ) - 1. (25) 

We define the fcth harmonic distortion by 

D h = lOlogdk, (26) 

and the total harmonic distortion by 

D = 101ogd. (27) 

III. CLIPPING MODEL 

Assume that the distortion may be represented as abrupt clipping of 
the positive peaks. Then 

v{d) = w cos (0,t) - w (cos (0,t)>„, (28) 

where w and r are parameters to be determined, and 

,„ , r O^0<t 

(0,t) = . (29) 

T<0<7T 

A simple calculation gives 

(cos (0,t))o = (t cos t — sin t)/tt. (30) 

The clipping factor is 

* = -20 log (cos t). (31) 

The Fourier coefficients are now obtained without numerical integration 
from the relations 

v\ = (w/tt){tt — t + - sin 2t) 

= w + . . . (32) 

i>k = [2w/ir(k 2 — 1)] (cos kr sin t — k~ l cos t sin kr) {k = 2,3, . . .) 

= -(2w/3t)t 3 + ...(k< t" 1 ), (33) 
and the total distortion from 

{v 2 ) B = (w 2 /2ir)(w + t cos 2t - - sin 2t) - m; 2 <cos (0,t)>I 

= (^/2)[1 + (4/15tt)t 5 + ...]. (34) 
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; = 



(37) 



For the determination of w,t consider (11) in the form 

v(B) = g-cos 6 + q+ - rj(u). (35) 

Comparing (28) with (35) shows that 

q- = w, q+ = -w (cos (d,T)) e , (36) 

and that j(v) is being approximated by 

r _1 (cos — cos t)<7_ 0< < t 

T<0<7T* 

The normalized waveform defined in (20) is being approximated by 

0(0) = cos (0,t). (38) 

When the model is valid, it is also valid to neglect j(a) in (7); it then 
follows from (7), (12) and (13) that w,t must obey the relations 

w = rg x (39) 

j(b)(l - cos r)" 1 = gl , (40) 

where 

b = w[cos t + (sin t — t cos t)/tt\. (41) 

Equation (40) requires that (37) be exact at = 0. 

The clipping model is suitable for use with a minicomputer since no 
numerical integrations are required. If/' (6) is easy to evaluate, the iter- 
ation of (40) is easy to do by trial and error. The accuracy of the model 
is best determined by comparison with the results calculated by the 
method of Section II. However, it is possible to obtain a corrected 
waveform in the region of clipping by solving the implicit relation 

j(»(B)) = (cos - cos r)g x (0 < 6 < r). (42) 

obtained from (37). 

IV. RECEIVER SENSITIVITY 

In Fig. 1 let the light power be 

U = U + U\ cos 

6 = wt, Uo^ux^- 0. (43) 

The photodiode assumed here is a specific n + xp + silicon unit having 
the dark current characteristic shown in Fig. 2. It is typical of a class of 
photodetectors developed by H. Melchior 7 for lightguide applications 
not requiring the sensitivity of avalanche photodiodes. 10 It has a quan- 
tum efficiency at 900 nm of r\ = 0.8 and a very low series resistance 
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Fig. 2 — Dark current-voltage characteristic of the photodiode. Parameters of the 
characteristic (44) for regions (t) and (it) are listed in Table I. 



« 1 fi. The characteristic of Fig. 2 has two exponential regions of the 
form 



;'(«/) = a(e0" - 1). 
The current generator g in Fig. 1 can be written 

g= yu. 



(44) 



(45) 



The parameters a,(5,y of the photodiode are listed in Table I. 

The earphone is the ring armature telephone receiver 8 which provides 
an essentially flat response over the voice band, 300-3300 Hz. We shall 
specify this earphone by a resistance p and power sensitivity 5, both as- 
sumed independent of frequency over the band. The response may be 



Table I — Photodiode parameters, 
t? = 0.80, 7 = eri/hi> = 0.58 V~ 1 (at 900 nm) 



(i) v < 0.37 V 



Hi) v > 0.37 V 



«(A) 
(a/3)" 1 (0) 



3.1 X 10" 12 
38.70 
8.3 X IO 9 



3.8 X IO" 9 
19.35 
1.36 X IO 7 
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characterized by the relation 

SPL = 81 + 10 log (p/8) (dB), (46) 

where SPL is the sound pressure level produced in a closed volume of 
6 cm 3 and p is the inband power available from a matched generator of 
resistance p. Equation (46) shows that, when p equals the power sensi- 
tivity 8, the telephone receiver produces a speech pressure level (SPL) 
of 81 dB, which is the average level found in the telephone network by 
surveys. 11,12,13 Representative values for p,d may be taken as 

p = 128 Q 

8 = 2.5 X 10~ 7 W. (47) 

The transformer is assumed to be ideal over the band with primary 
taps to give a variable turns ratio n. The dc resistance is assumed neg- 
ligible; this reduces the dc bias on the junction to zero, the most advan- 
tageous value. The inductance of the secondary must exceed 0.07 H 
determined by p and the low-frequency cutoff. It follows that the size, 
weight, and cost of the transformer would be approximately proportional 
to the maximum required value of n. In the following we specify the re- 
ceiver by n; the effective ac load resistance of the photodiode is 

r = n 2 p. (48) 

The sensitivity of a receiver of given n will be defined as the value of 
u\ which produces SPL = 81 dB at the ear; thus the sensitivity s is de- 
fined by the relations 

8 - i*i, pi = 8 (49) 

where pi is defined by (23). There is a range of optical powers u\ (and 
hence a range of sensitivities s) allowed by (49) because an increasing 
turns ratio n can be used to compensate (until nonlinearities become 
significant) for a decreasing optical power u\. Figure 3 shows s plotted 
versus n. At the minimum we find by the method of Section II the 
values 

n = 48, r = 0.29 Mfi (s = X ) (50) 

and 

X s «min = 2.6/lW. (51) 

The clipping model of Section III gives the same asymptotic straight line 
and the minimum value 2.5 p W. 

The asymptotic straight line represents the sensitivity s of the 
distortion-free receiver that would result if we could take j(u) = 0; the 
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Fig. 3 — Sensitivity s defined in (49) versus turns ratio n of transformer. The straight 
line shows the distortion-free receiver (52). 



equation of the line is 



y 2 rs 2 y 2 n 2 p A . 

8 = = s^ 

2 2 



(distortion-free receiver). 



(52) 



V. RECEIVER HARMONIC DISTORTION 

We define the amplitude level U\ by 

Ui = 20 log (ui/x). 



(53) 



The need for the factor of 20 in this equation in comparison to the factor 
of 10 in (46) results from (i) the photoelectric effect in the photodiode 
by which the electrical signal power produced is proportional to the 
square of the optical power modulation u\ and (ii) the desire to express 
U\ on the same logarithmic scale as SPL. In analogy to the present me- 
tallic telephone network we assume that the level reaching any receiver 
can be regarded as a random variable with a distribution 1 1>12 - 13 that is 
approximately normal with a variance of 7.8 dB. We define the dynamic 
range T of the optical channel 

r-C/ -<C7i>N. (54) 
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Fig. 4 — Total harmonic distortion D defined in (27) versus input signal level U\ defined 
in (53) relative to sensitivity level S defined in (56) for two values of n. Various sound 
pressure levels SPL are indicated; for all n,U\ = S corresponds to SPL = 81 dB. Points (t) 
and (it) are chosen for waveform examination in Fig. 5. 



where < U\)n is the average over the level distribution N (see Appendix 
A), and 



U u 20 log (u /x) 



(55) 



is the clipping level of the channel. This will always be a bottom-side 
clipping level; we will also assume for simplicity that it is a top-side 
clipping level. We assume T is a system constant maintained by the 
electrical -optical interfaces, whereas U and (Ui)n fall off with the 
transmission distance x of the particular lightguide. The dynamic range 
determines the probability that clipping will not occur, which we here 
call the quality. A more general discussion of T and speech quality is 
given in Appendix A. 

Harmonic distortion D&Ds, . . . and total harmonic distortion D are 
defined in (26) and (27) respectively. These quantities are functions of 
Ui and are only meaningful up to U\ = U = T + < Ui)n. We define a 
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Fig. 5 — Normalized waveform (20) for points (i) and (ii) of Fig. 4 for half of fundamental 
range of 0. The clipping model approximation is shown dashed and the parameter t is in- 
dicated on (ii). 

sensitivity level 

S = 20\og(s/ X )>0 (56) 

with a similar definition for § of the distortion-free receiver (52). For 
our calculations we have chosen to consider the receiver at the reference 
point x defined by 

(U 1 (x)) N = S, (57) 

or the distortion-free receiver at the reference point x~df defined by 

(U 1 (x DF )) N = S. (58) 

Since we will presently conclude that significant distortion in the receiver 
will be avoided in practice, the distinction between x and xdf need not 
concern us further. For illustrative purposes we choose the value T = 18.2 
dB. The reasonableness of this value in terms of speech quality is dis- 
cussed in Appendix A. Figure 4 shows D versus U\ — S for two receivers 
n = 47.7 and n = 13.4 out to U\ — S = 18.2. The steep rise of these curves 
and the absence of any extensive straight portions show that D is not 
dominated by the second harmonic but involves a large number of har- 
monics. At several points the SPL at the fundamental is indicated to show 
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Fig. 6 — Total harmonic distortion D (S + 18.2) at the channel clipping level for the case 
T = 18.2 dB versus sensitivity level S defined in (56). Dashed curve shows D(C) at receiver 
clipping level C defined in (61). 



that an 18-dB rise in level does not produce a corresponding rise in SPL 
when D is high. For a distortion-free receiver, (46) becomes with the help 
of (52), (53), and (56) 

SPE = 81 + Ui - S (distortion-free receiver). (59) 

The normalized waveform (20) is shown in Fig. 5 for two cases iden- 
tified as points (i), (ii) in Fig. 4. The distortion of the waveform is shown 
to be a gradual clipping of the positive peak. In the clipping model this 
is approximated by abrupt clipping out to 9 = r as indicated for curve 
(ii). We have obtained good results in calculating D%, D3 and D by the 
clipping model when D > —40 dB; higher harmonics are given with di- 
minishing accuracy. 

Figure 6 shows the total distortion at the channel clipping level 
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Fig. 7 — Clipping factor (21) versus level (53) for two illustrative values of n. Dashed lines 
represent abrupt clipping approximation (61) which defines the clipping level C of the 
receiver. The clipping discount A (see Appendix A) is chosen as A = 6 dB. 

D(S + 18.2) versus S for a receiver at the reference point. Conceivably 
an optimization could be based on an upper-limit objective for D at the 
channel clipping level. The dashed curve shows D(C) at the receiver 
clipping level (to be defined in the next section). Notice that D(C) is 
approximately constant, so in the present instance the quality matching 
principle is in effect equivalent to requiring D < —32 dB at the clipping 
level. 



VI. RECEIVER CLIPPING LEVEL 

The clipping factor has been defined in (21) for the general analysis 
and in (31) for the clipping model. The validity of the clipping model has 
been confirmed from the waveform and from calculations of D. Using 
the clipping model, we have calculated $(Ui) for various receivers n. (For 
brevity, n is not explicitly indicated in writing <I>.) Figure 7 shows the 
results for n = 13.4 and n = 6.5. It is known 14 that telephone speech 
quality, as determined in subjective listening tests, is not degraded ap- 
preciably by small clipping, <£ < A, where we call A the clipping discount 
and adopt the value 

A = 6 dB. (60) 

(This value has been deduced by us from an examination of the un- 
published work of A. M. Noll. 14 ) The discount is shown as a dashed line 
in Fig. 7. At each intersection of the discount with * we draw a line of 
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unit slope as shown. This represents abrupt clipping 

[0 Ui < C 



*(C/i) 



Ui - c u x >c 



at the receiver clipping level C. Thus C (as a function of n) is defined 
by the relation 

$(C + A) = A. (62) 

Notice that C is not very sensitive to the value chosen for A; any value 
of A in the range 2 to 8 dB would give about the same C(± 1 dB). 

If the concept of a clipping level is valid, the actual $(Ui) for the re- 
ceiver can be replaced with the abrupt clipping approximation (61). A 
test of the validity of the concept is the distortion D at U\ = C; for abrupt 
clipping D would be zero up to U\ = C. Figure 6 shows D(C) versus S; 
it is approximately constant around —35 dB. The departure of $(£A) 
from abrupt clipping above A is not of great significance, because the 
important question in quality determination is whether degradation 
occurs, not how much degradation has occurred. 

VII. QUALITY MATCHING 

The optimum sensitivity is the smallest value consistent with the 
requirement that the quality of the channel not be degraded by the re- 
ceiver. This gives the principle of quality matching expressed by the 
relation 

C = U (quality matching), (63) 

that is, the equality of receiver and channel clipping levels. This matching 
is to hold at all points x in the optical loop system. 

Figure 8 shows a schematic diagram of levels in a system versus the 
distance x from the optical source at the electrical-optical interface. We 
define the power transmission loss of the lightguide in the usual way 

TL(x) = 10 log [u(0)/u(x)] (dB). (64) 

At x = the interface injects the optical power 

u(0) = h + hicos 9 (65) 



at the levels 



It follows that 



tf = 201og(/i o /x) = t/ o (0) 

H 1 = 201og(/i 1 /x) = t/ 1 (0). (66) 



U (x)= Ho -2tl(x) 

Ui(x)=H 1 -2tl(x). (67) 
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DISTANCE x 



Fig. 8 — Levels versus distance in a quality matched lightguide loop system (receiver 
portion only). The transmit levels (referring to the optical source) are the channel clipping 
level U , mean amplitude level < U\ ) N , and transmission loss TL defined in (55), (54), and 
(64) respectively. Receiver sensitivity level 5 is also shown on the transmit scale. Receive 
level (68) is also shown. Reference point x is defined in (57). The optical source is charac- 
terized by H ; {Hi)n is defined in (66). 



In Fig. 8, TL, Uq, and (Ui)n are called transmit levels and are referred 
to the scale on the left. Also shown referred to the left scale is the receiver 
sensitivity level S. The receive level Rh(x) defined by 



RL(x) a (SPL(x)) N 



(68) 



is shown referred to the scale on the right. The reference point defined 
by (57) is denoted by x. The dynamic range T defined by (54) is a con- 
stant of the system. 

From (54), (57), and (63) we obtain the quality matching relations 



where 



C(x) - S(x) = r 

C(x) - C(x) = 2A(x), 

A(x) = tl(x) -tlCc). 



(69) 
(70) 

(71) 



Figure 9 shows C, C — S, and 5 plotted versus n. The implementation 
of quality matching is illustrated for the case r = 18.2 dB. We find point 
(i) on the C — S curve according to (69), which determines point (ii) on 
§ and (Hi) on C. The optical power at the source is then determined by 
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Fig. 9 — Receiver clipping level C defined in (61) versus n. Also shown are C — S,S, and 
S, where S is the sensitivity level (56) and $ is the sensitivity level for the distortion-free 
receiver defined by (52). Points (i) to (i>) trace the quality matching procedure based on 
(54), (63), (69), (70), and (74). 



(Hi) through the relation 



H = C(x) + 2TL 



(72) 



in terms of TL = TL(x). The choice of TL determines the average RL of 
the system, which is a system objective (not necessarily 81 dB) we must 
leave open at this time. We find point (iv) on C according to (70), which 
determines the optimum value of n for a receiver at x. It also determines 
point (v) on S. (The intersection with S, when it differs from §, is not 
of interest.) From (22) and (v) we determine the quantity 

P(x)=S(x)-S(x). (73) 

From (59), (67), and (68) the receive level for a distortion-free receiver 
is 



RL(x) = 81 + P(x) - 2A(x). 



(74) 



This is an excellent approximation for RL, because quality matching 
guarantees that distortion is too small to have any effect on response. 
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In Fig. 9 the lines C and § are given by the relations 

C = 70.5 - 44.3 log n (75) 

§ = 32.6 - 20 log n. (76) 

If T > 5 dB, it is justified to approximate 

C -S = 37.9 - 24.3 log n. (77) 

It follows that 



C(x) = 1.4 + 1.82 T 


(78) 


P(x) = 0.90 A(x) 


(79) 


log n = 1.56 - 0.041 T + 0.045 A. 


(80) 


Thus the receive level is 




RL « Rt = 81 - 1.1 A. 


(81) 


and the clipping level at the optical source is 




H = 1.4 + 1.82 T + 2TL. 


(82) 



VIII. SUMMARY AND DISCUSSION 

We have obtained a practical algorithm suitable for a large computer 
for solving the nonlinear integral equation (3) referring to Fig. 1. The 
equation implies an integral restraint (4) on the solution which is an 
unusual feature that removes it from the types found discussed in texts 
on nonlinear networks. The load r is taken as zero at dc and a constant 
resistance at all signal frequencies of interest. This implies that the load 
is dispersive at frequencies below a certain cutoff frequency (300 Hz in 
the receiver problem). Ordinarily a nonlinear dispersive circuit requires 
nonlinear differential equations to describe it. Here we have avoided the 
differential equation and obtained instead an integral equation (3) by: 
(i) asking only for the periodic solution, and (ii) treating separately the 
ac and dc voltages with different values of r. The assumption of zero dc 
resistance is convenient and usually appropriate in practice, but the 
analysis presented in Section II can easily be generalized to any value 
of dc resistance. The nonlinear conductor j(») is passive (/(0) = 0) and 
monotonically increasing (/(")' > 0) but otherwise arbitrary. The con- 
ditions on j{v) rule out negative-resistance instabilities and guarantee 
a unique solution. 

The method presented in Section II is exact in principle; it is based 
on the fact that any functional of the solution (e.g., a Fourier coefficient) 
can be calculated explicitly by integration once two parameters (a, b) 
have been determined from the implicit relations (12), (13). Standard 
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routines are available for solving simultaneous implicit relations to any 
desired precision. In using this method we have usually obtained satis- 
factory convergence with no special precautions. When convergence 
problems are encountered, the answer is to start the iteration with better 
estimates for a,b. 

The clipping model described in Section III was originally worked out 
to provide initial estimates of a, b for use in the exact method. It soon 
became apparent, however, that it is sufficiently accurate in the receiver 
problem for all calculations of sensitivity, total distortion (when D > -40 
dB), and clipping factor. The reason for this is that the exact waveform 
comes close to the clipped waveform assumed in the model. The clipping 
model is not recommended for the calculation of specific harmonics 
higher than the second. Generally the clipping model is expected to be 
useful whenever j(v) is a rapidly increasing function. In this model dis- 
tortion (clipping) is represented by a single parameter t which is de- 
termined from the implicit relation (40). No numerical integrations are 
involved in using the model, which makes it convenient for use with a 
minicomputer. 

The analysis of the photovoltaic receiver in the remainder of the paper 
is based on a sinusoidal input waveform. The calculation of the sensi- 
tivity s, defined in (49) as a measure of response based upon producing 
a certain reference sound level at the fundamental of 81 dB, is presented 
in Section IV. The photodiode and the earphone assumed in this cal- 
culation are the best presently available for the purpose. The transformer 
secondary must have an inductance of at least 70 mH, so the turns ratio 
n is adjusted by means of taps on the primary. By adjusting n, any value 
of s down to the minimum 2.6 juW can be obtained. However, the size and 
cost of the transformer are expected to increase approximately as the 
maximum value of n. 

The minimum s = x shown in Fig. 3 is a nonlinear effect having 
nothing to do with the impedance matching concept of linear circuit 
theory. The small signal resistance of the junction is given in Table I, 
(ajS)" 1 = 8.3 X 10 9 G which is larger than r(x) = 2.9 X 10 5 fi by a factor 
« 3 X 10 4 ! This shows the necessity for a nonlinear analysis. For n > 48 
the distortion increases very rapidly at the reference level assumed for 
the calculation of s. This does not mean, however, that n > 48 is ruled 
out for receivers operating at much lower levels. In a properly designed 
system the receiver must respond almost like the distortion-free receiver 
defined in (52) at the levels to which it is subjected. The reference level 
SPL = 81 dB is the overall average level for receivers in the existing 
telephone system. 

At a particular point, the reference point, in a loop system the receive 
level defined in (68) will equal 81 dB, and at other points its value will 
depend on the transmission loss from the reference point to that point. 
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For a discussion of total harmonic distortion D we have considered a 
receiver at the reference point in Section V. We find (Fig. 4) that D rises 
very rapidly with level up to about D « — 20 dB and then begins to bend 
over. The value of D of interest is the value at the channel clipping level 
shown versus sensitivity in Fig. 6 for an assumed dynamic range of 18.2 
dB. The results confirm that D is much too high in the minimum sensi- 
tivity receiver for use at the reference point. It is conceivable that some 
criterion on D (e.g., D < —20 dB) could be used as the basis for optimizing 
the receiver (choosing n). However, we believe that clipping provides 
a more objective and less arbitrary basis for optimization. The existence 
of clipping is shown by the waveform of Fig. 5 and by the good agreement 
between the clipping model and the exact method of D. 

The clipping factor defined in (21) has been calculated as a function 
of level in Section VI. At a certain value A, called here the clipping dis- 
count, abrupt clipping begins to degrade speech quality. 14 Therefore an 
abrupt clipping approximation has been fitted to the receiver clipping 
factor at the value A. This defines the clipping level C as illustrated in 
Fig. 7 for the choice A = 6 dB. Actually C is not very sensitive to A. The 
total distortion at level C would be zero for abrupt clipping. We find D(C) 
« —35 dB for all receivers (Fig. 6); in our opinion this is small enough to 
confirm the validity of the clipping level concept. 

Quality matching as an optimizing principle is defined in Section VII. 
It amounts to setting the channel and receiver clipping levels equal as 
in (63). The quality of speech transmitted in an analog channel in which 
the only nonlinearity is abrupt clipping is determined by the dynamic 
range, defined for sinusoidal signals in (54). The dynamic range T is in- 
dependent of lightguide transmission loss so we consider it a constant 
of the optical loop system. In Fig. 8 we pass from a strictly device view- 
point to a system viewpoint. The quality matching relations (69) de- 
termine the optical power needed at the reference point as well as the 
optimized receiver and receive level throughout the optical loop system. 
The various levels are shown schematically in Fig. 8 as a function of 
distance x from the optical source assuming the same source power for 
all loops. The curves required to obtain the solution are shown in Fig. 
9, and approximate equations for the solution are given in (78) through 
(82). 

In the existing system the mean loop insertion loss 15 is about 5 dB. The 
receive level defined in (68) is 81 dB at a point 5 dB from the central 
office, and 81 dB is approximately the mean loop receive level. This is 
not necessarily an objective for future loop planning, so we here em- 
phasize that our analysis contains no such assumption. The use of 81 dB 
as a reference is only a convenience and involves no loss of generality. 
The reference point x in Fig. 8 with the transmission loss TL need not 
be the "average point" for the loop system. The mean loop receive level 
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is from (81) 

<RL> L = 81 + 1.1 (TL - (TL> L ) (83) 

where ( )i denotes an average over loops. This can be set at any desired 
level by properly choosing TL. 

The quality matching concept requires a receiver with adjustable turns 
ratio which probably involves some added cost compared with a fixed 
receiver. However, any fixed receiver would give a receive level varying 
as-2A, 

RL = 81 - 2 A (fixed receiver), (84) 

whereas quality matching gives (81) varying approximately as —A. This 
effect is shown in Fig. 8 by the different slopes of RL and (U\)n. Clearly 
this is a desirable effect which utilizes the capabilities of the receiver to 
the fullest extent and permits serving a radius approximately twice that 
which would be possible with a fixed-sensitivity receiver having the 
sensitivity of an optimized receiver at the reference point. The actual 
range of the loop system would probably be limited by objectives on the 
minimum RL which we are not discussing here. Another limitation which 
we can only mention is that of transformer cost. From (80) we find that 
the cost, as measured by n, doubles for an increase of 6.7 dB in TL. 

To illustrate the principles being discussed, we have chosen a realistic 
case specified by 

T = 18.2 dB 

TL = 5 dB (illustrative) 

A(x) = 7.5dB. (85) 

Points (£) to (u) in Fig. 9 trace the solution for this case. At the reference 
point x we find 

C(x) = t/ (Jc) = 34.6dB 

u o 0t) = 0.14mW 

s(x) = 17mW. (86) 

At the point x we find 

C(x) = 19.6 dB 

n(x) = 14.1 

RL = 72.8 dB. (87) 
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At the source we find 

H = 44.6 dB 
(H 1 ) N = 26.4 dB 
h = 0.45 mW 
peak source power = 2 h = 0.9 mW. (88) 

The fact that the peak source power is 2 h follows from our assumptions 
of equality of top and bottom side clipping and of occurrence of bottom 
side clipping at h = 0. 

The sensitivity at the reference level s(x) = 17 nW is 50 times smaller 
(better) than that of the best optoacoustic receiver, the xenon photo- 
phone 4 (s = 0.9 mW). Furthermore, the photophone is a fixed receiver. 
Therefore the photovoltaic receiver is clearly superior for loop applica- 
tions. 

The "system" referred to here should be regarded as a relatively small 
subsystem of the loop plant serving a special class of customers who have 
lightguides running to their premises primarily to provide high capacity 
services. At the central office or other junction point there must be an 
electrical-optical interface containing the optical source for the receiver. 
The peak power of this source according to (88) should be 0.9 mW. This 
power into the lightguide is within the capabilities of present day het- 
erojunction laser diodes 16 but about an order of magnitude above the 
capabilities of luminescent diodes. 16 To a limited extent the power at 
the source might be varied to compensate for the transmission loss; one 
could imagine all adjustment being done at the source instead of the 
receiver. Although we have described a system with fixed source power 
in Fig. 8 and believe that this is the most likely type of system, no change 
is required in the equations to treat the variable source. The dynamic 
range, however, would be determined by the interface circuitry and held 
to a uniform value to maintain transmission quality. 

It may be objected that our nonlinear analysis has been based on a 
sinusoidal input signal whereas a telephone is required to transmit 
speech. Appendix A contains a discussion of the extension of the theory 
to a general waveform and the general definition of dynamic range. It 
is argued that the clipping level is valid for any waveform. A quality Q 
is objectively defined on the basis of the approximately normal distri- 
bution of speech levels in the telephone system. Figure 10 shows Q versus 
T for speech and for a hypothetical system having sinusoidal signals. 
Essentially, Q is a probability that speech is received without degradation 
caused by clipping. The illustrative case T = 18.2 dB used in this paper 
corresponds to Q = 86 percent and Q coa = 99 percent. An objectively 
defined quality is not equivalent to a grade of service determined sub- 
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Fig. 10 — Objective speech quality Q defined in (104) versus dynamic range T defined 
in (93) plotted on probability scale. The corresponding quantity for sinusoidal waveform 
is shown dashed. 

jectively from listening tests, although we would expect a strong corre- 
lation between the two. 

We shall not speculate on the circuitry of the interface, but we offer 
an interesting observation in Appendix B. On the basis of the light- 
emitting characteristic of laser diodes and the signals from the metallic 
network existing at the central office, it is shown that the direct appli- 
cation of the metallic network signal current (no amplifiers or trans- 
formers) to the laser would produce a somewhat greater level than (Hi)n 
called for in (88). 

Let us suppose that TL = 5 dB in the illustrative case (85) is the av- 
erage loss and the lightguide attenuation constant 17 is 2 dB/km. Then 
the reference point, which is also the average point, is at the distance 



x = 2.5 km = 8.3 kft (illustrative). 



(89) 



This distance may be compared to the average loop length 15 10.3 kft in 
the present metallic loop plant. A more conservative estimate of atten- 
uation, like 6 dB/km, reduces (89) by a factor of three. Thus, the optical 
loop system seems to be limited to a somewhat smaller radius of service 
than the metallic loop system. 
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APPENDIX A 

Speech Quality and Dynamic Range 

We continue to let 6 denote a normalized time variable, but in (43) and 
(65) we replace cos with a more general normalized waveform \p(6) 

max tffi) = -min i//(0) = 1 

W)>* = 0. (90) 

Define the peak factor 

^ = -201og^ rm5 (91) 

and signal level 

U= 20 log (un/w/x) 

= U l - * (92) 

where U\ is given by (53). Define the dynamic range 

r=(201og(u o /2 1 / 2 ui^ rms )) N 
= U -(U) N -S (93) 

where the channel clipping level U is given by (55) and where 

(U) N = J^mWdN, (94) 

N being a distribution function. For \p(6) = cos we have ^ = 3 dB and 
(93) reduces to the definition first given in (54). For the distortion-free 
receiver the SPL is 

SPi; = 81 + t/-.S+3, (95) 

which reduces to (59) when ^ = 3. We assume (61) and (63) remain valid 
but write (61) in terms of U 

K } u+*-u u>u -*. 

The levels U are distributed 11,12 - 13 such that the probability that 
U<Xis 

P{U <X) = N[(X - X )/ff] (97) 
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where 

X = (U) N , <r=7.8dB (98) 

and N(x) is the normal distribution 

N(x) = (27T)" 1 / 2 C* e-W dt. (99) 

It is known from subjective studies 14 of the effect of clipping on telephone 
speech quality that no degradation occurs for $ < A, where we take 

A = 6 dB. (100) 

(This is a deduction from the work of A. M. Noll for which we assume 
responsibility.) The probability that $ < A is 

P(* < A) = N[(U + A - * - (U) N )/a) 

= N[(T + A + 3 - ¥)/4 (101) 

The speech waveform is characterized by the value 18 

tf = 18.6 dB (speech), (102) 

so that 

A + 3 - * = -9.6 dB (speech). (103) 

The objective definition of quality is 

Q = 100 X P($ < A) (speech) 

= 100 X N[(T - 9.6)/7.8]. (104) 

A corresponding definition can be given for sinusoidal signals having the 
same level distribution assuming A = and ^ = 3 

Q coa = 100 X iV(r/7.8). (105) 

Figure 10 shows Q and Q cos versus T on a probability scale. 

APPENDIX B 

Comments on the Optical Source 

A laser diode has a threshold current for lasing and at higher currents 
a steeply rising emission h as a function of current i. In the lasing region 
we assume 19 

dh/di = e 

e = 2.3 W/A. (106) 

If a sinusoidal current of amplitude ii flows through the laser (super- 
posed on a suitable bias current), the optical signal level Hi defined in 
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(66) is 

ffi-201og(s'i/x). (107) 

Suppose that the current i\ is that supplied to a matched load by a 
generator of resistance r and available power p; then (107) can be writ- 
ten 

H r = 20 log [2 1 ' 2 ( € /x)(p/r) 1/2 ]. (108) 

Finally, suppose that the generator is the central office and p corre- 
sponds to the mean signal level; representative values are 20 

p = 2/tW, r = 1166 ft. (109) 

It then follows that the mean optical amplitude level at the source is 

<>7 1 >7v = 34dB {hi = 0.13 mW). (110) 

This compares very favorably with the level called for, 26 dB, in (88). 
This shows that the current flowing in a matching resistor (1166 ft) 
combined with the "gain" of the laser emission characteristic provides 
more than enough signal without a matching transformer in the cir- 
cuit. 

If a matching transformer is used, r in (108) is replaced by the dynamic 
laser diode resistance 19,21 /•<* » 1.5 ft, giving hi ~ 3.8 mW. For a lumi- 
nescent diode an appropriate value of e in (106) is 16 c « 5 X 10~ 4 W/A 
giving for the transformer matched case with r^ = 1.5 ft the power hi « 
0.8mW,(H 1 ) w = -10 dB. 
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